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Abstract. 

The effects of small-scale heterogeneity in land surface 
characteristics on the large-scale fluxes of water and energy in the 
land-atmosphere system has become a central focus of many of the climatology 
research experiments. The acquisition of high resolution land surface data 
through remote sensing and intensive land— climatology field experiments (like 
HAPEX and FIFE) has provided data to investigate the interactions between 
microscale land-atmosphere interactions and macroscale models. One essential 
research question is how to account for the small scale heterogeneities and 
whether ’effective’ parameters can be used in the macroscale models. To 
address this question of scaling, the probability distribution for 
evaporation is derived which illustrates the conditions for which scaling 
should work. An correction algorithm that may appropriate for athe land 
parameterization of a GCM is derived using a 2nd order linearization scheme. 
The performance of the algorithm is evaluated. 



Introduction. 


Research into land-atmospheric interactions suggest a strong 
coupling between land surface hydrologic processes and climate (Charney et 
al., 1977; Walker and Rowntree, 1977; Shukla and Nintz, 1982; and Sud et 
al.,1990.) Due to this coupling, the issue of 'scale interaction' for land 
surface-atmospheric processes has emerged as one of the critical unresolved 
problems for the parameterization of climate models. To help resolve this 
issue, the understanding of the scaling properties of water and energy fluxes 
with their corresponding storage terms (especially soil moisture) has been an 
important scientific objective of land-climatology experiments like FIFE (see 
Sellers etal. 1988) and GCIP (see WCRP, 1992). 

In fact, the acquisition of high resolution land surface data 
through remote sensing and intensive land— climatology field experiments (like 
HAPEX and FIFE) has provided data to investigate the interactions between 
microscale land-atmosphere interactions and macroscale models. One essential 
research question is how to account for the small scale heterogeneities and 
whether 'effective' parameters can be used in the macroscale models. The 
current scientific thinking on this issue is mixed. For example Sellers et 
al. (1992) claim that analysis of the FIFE data supports that 
land-atmospheric models are almost scale invariant, a conclusion also reached 
by Noilhan (personal communication) using HAPEX— MOB I LHY data. Counter 
arguments have been made by Avissar and Pielke (1989) who found that 
heterogeneity in land characteristics resulted in sea-breeze like 
circulations, and significant differences in surface temperatures and energy 
fluxes across the patches. The results of analyses presented later in this 
paper suggest that soil moisture is the critical variable that controls the 



non-linear behaviour of land-atmospheric interactions, and the effect is most 
pronounced when soil moisture heterogeneity is such that part of the domain 
is under soil-vegetation control and part is under atmospheric control. 

The understanding of the scale problem is critical for new climate 
projects such as the Global Energy and Water Experiment (GEWEX) or the Earth 
Observing System (Eos). It is also important for the parameterization of the 
■microscale land-surface hydrology in climate models, and crucial in our 
understanding in how to represent sub-grid variability in such macroscale 
mode 1 s . 


From a modeling perspective, it’s important to establish the 
relationship between spatial variability in the inputs and model parameters, 
the scale being modeled and the proper representation of the hydrologic 
processes at that scale. Figure 1 presents a schematic for modeling over a 
range of scales. Let us consider this figure in light of the terrestrial 
water balance, which for a control volume may be written as- 


<f§.> _ <P> _ <E> - <Q> (1) 

ot 

where S represents the moisture in the soil column. E evaporation from the 
land surface into the atmosphere. P the precipitation from the atmosphere to 
the land surface, and Q the net runoff from the control volume. The spatial 
average for the control volume is noted by (•). 

Equation (1) is valid over all scales and only through the 
parameterization of individual terms does the water balance equation become a 
‘distributed’ or ’lumped’ model. By ‘distributed* model, we mean a model 
which accounts for spatial variability in inputs, processes or parameters. 



Lumped Versus Distributed Models. 

Figure 1 presented a framework for considering the relationship 
between distributed and lumped models. Wood et al. (1988). Wood et al. 
(1990), Wood and Lakshmi (1993) studied the behaviour of aggregated inputs 
and the resulting hydrologic responses which lead to the concept of the 
representative elementary area, a scale where a statistical representation 
can replace actual patterns of variability. In this paper we compare the 
output between a macroscale, distributed model and a lumped model to try and 
determine when the macroscale model provides an accurate response when 
compared to the average of the distributed model. 

Two distributed model are used in latter sections: one is based on 
the model described in Famiglietti and Wood (1992) (referred to as TOPLATS - 
Topgraphic-Land-Atmospheric-Transfer-Scheme) and is an extension of the model 
described in Famiglietti et al . (1992); these models have been applied to the 
Intensive field campaign periods (IFCs) during FIFE of 1987. TOPLATS can 
include variability in topography, soils, net radiation and vegetation. The 
first two, topography and soils, leads to variations in soil moisture under 
the TOPMODEL framework; the latter two lead to variations in potential and 
actual transpiration. The second model (Wood et al., 1992) described here 
has variations in infiltration capacities across a catchment of GCM grid 
square and is referred to as the VIC model. This model has been used in GCM 
climate simulations (see Stamm et al , 1993) and is the model used in this 
paper to demenstrate the impact of a second order correction term to the 
lumped model to account for sub-grid heterogeniety . 

Derived Distribution of Soil Moisture and Evaporation. 

The concept of the representative elementary area leads to a 
statistical description of the sub-grid variability in water table depths. 


soil moisture, evapo transpiration, and so forth. Models at this scale are 
referred to as macroscale, distributed models. For TOPLATS, the distribution 
in the soil-topographic index leads to a distribution in water table depths 
(see Beven and Kirkby, 1979; Wood et al., 1990; or Famigliettl et al . . 1992). 
From the distribution of water table depths z. and the soil characteristic 
relations which relates the soil matrix head. to soil moisture 0 as a 
function of soil properties, the statistical distribution for surface soil 
moisture can be derived. Certain assumptions are usually applied, for 
example a steady state vertical flow which leads to cty/3z = 1. Actual 
evaporation e & (or transpiration from vegetation) depends on the availability 
of soil moisture (i.e. a soil controlled rate) and the atmospheric demand for 
moisture; the actual evaporation being the minimum of the two rates. Thus 
given the atmospheric demand and the the statistical distribution of soil 
moisture, the distribution of the actual evaporation can be derived. For 
simple functional forms the mapping of z -* 0 -* e can be done analytically; 
in any case it can be done through simulation. Figure 2 provides some 
results using TOPLATS for two conditions' Figure 2a is for quite dry 
conditions - low water table - for two times during the day. Figure 2b is 
for the same dirurnal times but for a wetter (but not extremely wet) 
condition. The parameters for the curves are taken from Famigliettl and Wood 
(1992) and represent conditions for Kings Creek area of FIFE in Kansas. The 
figures have been divided into four panels that show the derivation of the 
bare soil evaporation distribution. Panel (i) in the lower right corner 
gives the probability distribution for water table depths derived from the 
soil- topographic index of TOPLATS. and soilwater table depth to surface soil 
moisture as discussed ewater table depth to surface soil moisture as 
discussed earlier in this section. These two figures could gewater table 
depth to surface soil moisture as discussed earlier in this section. These 



two figures could generate a derived probability distribution for surface 
soil moisture. This is not shown here. The upper left portion of the figure 
(iii) gives the relationship between surface soil moisture and actual 
evaporation for the two times during the day. The maximum evaporation rate 
is the potential rate, which is lower during the early morning and late 
afternoon. For portions of the catchment where surface soil mositure is 
high, the actual evaporation rate is equal to the potential rate. For drier 
areas, the rate is lower. The resulting probability distribution for the 
actual evaporation is shown in the lower left portion of the figure. This 
panel is divided into two, the top giving the distribution for the time 
related to the lower potential evaporation rate. 

Inspection of the derived distributions for the two times and two 
conditions reveals that for the very dry conditions the distribution for 
evaporation is narrow and the average water table depth can be used to 
estimate the average evaporation rate. This is because the soil moisture - 
evaporation function is essentially linear in the range of soil moistures 
representing the dry conditions. For the wetter condition, the function is 
non-linear and the range of soil moistures contain areas which are very dry 
(having low evaporation rates) and wet (having rates at potential.) If the 
conditions were even wetter, then the distribution of evaporation rates would 
be at the potential rate. It is at these intermediate conditions where the 
non-linearities appear to have the greatest impact. 

Comparisons between the average bare soil evaporation and that 
estimated using the average depth to the water table are given in Table 1 for 
both conditions and the two times. From Figure 2 and Table 1 it is clear 
that dry conditions during periods of high atmospheric demands result in 



evaporation rates that are biased low. IXiring periods with high soil 
moisture (or extremely low soil moisture) or during early or late times 
during the day when the atmospheric demands are low, the bias from using the 
average water table depth is minimal. 

To test the sensitivity due to dry soil conditions and to compare 
the distributed water-energy balance model (TOPLATS) to a lumped 
representation (one-dimensional model or a first order model), comparisions 
were made between the models for 5 days during the October 1987 FIFE 
intensive field campaign, IFC-4. This period had the driest conditions 
observed during the 1987 experiment. Figure 3 shows the simulations for 
October 5-9, 1987. The models were run at a 0.5 hour time step to capture 
the diurnal cycle in potential evapotranspiration. Three models are 
compared: a fully distributed model, a macroscale model in which the spatial 
variability is considered statistically and a lumped one-dimensional model in 
which parameters and inputs are spatially constant. 

The one-dimensional model predicts well the evapotranspiration 
during the morning and late afternoon when the atmospheric demand is low, but 
fails to accurately predict this flux during the middle portion of the day 
when soil and vegetation controls limit the actual evapotranspiration. It is 
during this period that the sensitivity is high and by ignoring the spatial 
variability in soil moisture the lumped model serverely underestimates the 
catchment-scale evapotranspiration. During wet periods, the one-dimensional 
model nay work quite well. This complicates the linkage between a 
distributed and lumped representation since the appropriateness of the 
simpler representation varies with the state of the system. 



Linearized. 2nd Order Model for Sub-grid Variability. 

While these results imply that distributed models are needed to 
accurately account for sub-grid variability in soil moisture and the 
resulting evapo transpiration, such models may be computationally burdensome 
when incorporated within a GCM. An alternative approach would be to develop 
correction schemes for the often used * lumped models. Such a scheme should 
have correction terms that vary with soil moisture conditions. 


A lumped representation (or what will also be referred to as a 
one-dimensional representation) is obtained by using spatially constant 
values for 'state' variables; in the case of TOPLATS this would be the 


soil- topographic index and vegetation parameters; in the case of VIC the soil 
moisture. The effect of representing the distributed model by a lumped 
model, or equivalently by replacing the spatially variable parameters and 
Inputs by average values, will depend on nonlinearities in the model. 
Conceptually this can be seen by considering a second order Taylor's series 
expansion about the mean for the function y = g[x.0] where 0 are fixed 
parameters and x variable with mean p(x) and variance ct(x) . A first order 
approximation for y is Pj(y) ^ g[fi(x),0], while a second order approximation 
would be 


n 2 (y) - gO(xM] + \ 



o(x) 

M(x) 


(3) 


Differences between Pj(y) and p 2 (y) depend on the magnitude of the 
second term in equation (3) — the sensitivity term. As an illustrative 
example, consider the estimation of downslope subsurface flows, qj. within 
TOPLATS with and without considering variability in the local water table z^ . 
TOPLATS relates q A to Zj by q A = Tj tan/3 exp(-f Zj). Thus a first order 
approxination of the mean subsurface flow would be 



(4) 


= Tjtan/3 exp(-f z) 

while a second order approximation would be 

P 2 ( q i) = Tj tan/3 exp(-f z) + | { ^ tan/3 f } 2 exp(-f z) a( Zj) 

If we scale ^(qj) by ^l^i^ 01141 reco 8 nlze that 

2 aT e 

CT ( z j) - f a ( ln Tj tan/3 ^ 


(5) 

( 6 ) 


we obtain 


■VV 


= 1 + 0.5 o(ln 


aT 

T tanP J 


(7) 


Analysis of the soi 1- topographic index for Kings Creek yields a variance of 
3.25. This results in the first order estimate for qj of being biased low by 
approximately 65%. Since the subsurface flows and the local water table are 
related and since the local water table depth effects the surface soil 
moisture which subsequently determines the soil evaporation and infiltration 
rates, it’s clear that the lumped model may very well lead to significant 
biases in the water balance fluxes. 


For the more complex functions used for bare soil evaporation and 
transpiration, the sensitivities can be determined through simulation. For 
these functions the sensitivities will change with the state of the catchment 
(wet or dry). For example Figure 4 gives the vegetation transpiration and 
soil exfiltration capacities used to model the FIFE data (Famiglietti and 
Wood. 1992). Notice that at low and high soil moisture values the 
transpiration capacity function is essentially linear and the sensitivity 
would be low to soil moisture variations in these ranges. For volumetric 
moisture contents in the range 0.2 - 0.3, the sensitivity of the 
transpirtaion capacity function is high. As can be seen from Figure 4, 
sensitivity characteristics for soil exfiltration capacity would be high for 
soil moisture values greater than about 0.3. 


moisture, which for any fractional area with capacity i greater than i Q can 
be estimated as 0 = i <) Vl • where tj is the soil porosity. Using the results 
for a Beta distribution between [i Q . l m 3 results in a mean and variance for i 


conditional upon iM of 


i = *o + 


Var(i) = (i m -i o r 


B 


(1+B)^(2+B) 


( 11 ) 

( 12 ) 


A second order mean soil moisture can be estimated using (7) as 


in in 

0 = + o Var(i) 

T i J 


(13) 


Var(0) = -ZTZ — Var(i) 0 4 ) 

i 

To estimate bare soil evaporation we can use a Philip’s form of the 

- 1/2 

exfiltration capacity (Eagleson. 1978) which is of the form E = 0.5 S e t 

where S is a desorptivity term that can be written in the form 
e 

S = K(0-0 ) G/2 . Here K and C are parameters which depends on soil 
e v r J 

characteristics. Using the mean and variance of 0 will yield a second order 
model for the mean soil evaporation rate whose sensitivity will depend on the 
variance of 0 and the sensitivity term for the evaporation function. After 
some simple algebra, the sensitivity term d^d© 2 can be written as 


A = ( | | -1 ) ^—2 

do 2 2 2 (0 - e r r 

where E is evaluated at (1, and 0 f is the residual soil moisture. 


(15) 


To test the algorithm and compare it with the distributed VIC 


model, simulations were run varying the initial soil water capacity (W q ) and 



for different shape parameters for the distribution of infiltration 

capacities, (see equation 8). The simulations used potential evaporation 

data from IFC-4 of FIFE’87. The following initial parameters were used: a 

maximum infiltration capacity. I* = 30.5 cm and a shape parameter of B = 0.3. 

Using an initial wetness corresponding to i Q = 10 cm yield the results shown 

in Figure 7. The solid line represents the distributed VIC model in which 

evaporation is estimated using 100 slices of the soil moisture distribution. 

The dotted line is the solution using the average value of soil moisture for 

where i > i (i.e. a 1st order linear model) and the dashed line is the 1st 

order model plus the correction term (i.e. the 2nd order, linearized model.) 

For this case, the correction term is about 50% of the estimate using just 

the average soil moisture. If the conditions are wetter than Figure 7. i.e. 

have i = 15 cm. then we get the results shown in Figure 8. Here the lumped 
o 

model does very well. If conditions drier than Figure 7 prevail, i.e. i Q = 5 
cm (which is very dry) then the linearized model does very poorly: as shown 

in Figure 9. 

The results do depend on the value of B, the shape parameter. For 
example having B = 1.3 and i Q = 10 cm (the conditions of Figure 7) resulted 
in much better performance of the correction algorithm as can be seen in 
Figure 10. As shown in Stairnn et al. (1993). the distribution of B ranges 
globally from about 0.3 to 2.5 suggesting that the applicability of such a 
correction algorithm may be widespread. Furthermore, it can be determined 
before hand where the algorithm should work, and under what soil moisture 
conditions. This suggests that for those GCM grid squares with sufficient 
moisture or favorable infiltration capacity shape parameters, the simple 2nd 
order algorithm can be implemented. For condition too dry. the distributed 
model can be run for those particular time steps. This approach would lead 
to the most efficient and accurate computational effort. 



Conclusions 


The effect of subgrid variability in soil moisture on evaporation 
has been investigated with the aim of resolving whether effective (or 
average) values for soil moisture can replace the distribution found within a 
catchment or GCM grid square. It appears that there is a critical range of 
Intermediate values for which the subgrid variability has a significant 
impact of grid total evaporation (and transpiration). This arises from the 
non-linearity between soil moisture and evaporation within this critical 
range, and the essentially linear behavoiour outside this range. 

This lead to an initial attempt in developing a 2nd order, 
linearized, model for evaporation that could be incorporated with GCMs. 
Initial performance of this algorithm is encouraging with the correction term 
representing about 50% of the evaporation predicted based on only using the 
average soil moisture value. For extremely dry conditions the linearized 
model still under estimates evaporation which may result in using the fully 
distributed model in these conditions. 
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Table 1. Average Evaporation Rates 


Average based on the 
Variable Evaporation 

Average based on using 
average z 

Soil 

Condition 

Potential 

Low 

Evaporation (irnn/hr) 
High 

Potential 

Low 

Evaporation (nm/hr) 
High 

Dry 

Wet 

0.051 

0.11 

0.037 

0.088 

0.054 

0.12 

0.046 

0.097 





Figure 4- Vegetation transpiration capacity and bare soil evaporation 
capacities as functions of soil moisture (from Famiglietti and Wood, 1992). 

Figure 5: (a) Actual soil evaporation for five different times during a day, 

notice that the actual level is the minimum of the capacity (figure 4) and 

the potential, (b) The evaporation sensitivity term /d&^ for the same 

s 

times as in (5a). 

Figure 6: The vegetation transpiration sensitivity term . 

Figure 7: Comparisons of model derived latent heat estimates using the 

distributed VIC model for medium dry conditions, a 1st order linear model and 
a 2nd order linear model. Pamameters for the simulation are given in the 
text. 

Figure 8: Comparisons of model derived latent heat estimates using the 

distributed VIC model for wet conditions, a 1st order linear model and a 2nd 
order linear model. Pamameters for the simulation are given in the text. 

Figure 9: Comparisons of model derived latent heat estimates using the 

distributed VIC model for very dry conditions, a 1st order linear model and a 
2nd order linear model. Pamameters for the simulation are given in the text. 

Figure 10: Comparisons of model derived latent heat estimates using the 

distributed VIC model for the moisture conditions used in figure 7 and a 
modified soil capacity shape parameter. Also shown are a 1st order linear 
model and a 2nd order linear model. Pamameters for the simulation are given 


in the text. 
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